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Many dark energy (DE) models have been proposed, in recent years, to explain acceleration of 
the Universe expansion. It seems necessary to discriminate the various DE models in order to check 
the viability of each model. Statefinder diagnostic is a useful method which can differentiate various 
DE models. In this paper, we investigate the statefinder diagnosis parameters for the holographic 
dark energy (HDE) model in two cosmological setup. First, we study statefinder diagnosis for HDE 
in the context of flat Friedmann-Robertson-Walker (FRW) Universe in Einstein gravity. Then, we 
extend our study to the DGP braneworld framework. As system’s IR cutoff we chose the Hubble 
radius and the Granda-Oliveros cutoff inspired by Ricci scalar curvature. We plot the evolution of 
statefinder parameres (r, s} in terms of the redshift parameter 2. We also compare the results with 
those obtained for statefinder diagnosis parameters of other DE models, in particular with ACDM 
model. 


I. INTRODUCTION 

Recent astronomical observations such as type la supernovae (SNIa) J|, large scale structure (LSS) Q, and the 
cosmic microwave background (CMB) anisotropies Q confirm that our Universe is currently undergoing a phase of 
accelerated expansion. DE which is responsible for this acceleration has anti-gravity nature and hence push the 
Universe to accelerate. Although the nature of such a DE is still unknown, many candidates have been proposed 
for DE in the literatures. The first and simplest candidate for DE is the cosmological constant with equation of 
state (EoS) parameter ui = — 1. But unfortunately, it suffers from some serious problems such as fine-tuning and 
coincidence problems. A comprehensive but incomplete list of DE scenarios include scalar fields (such as quintessence 
[ll|, phantom (ghost) field [•}•. K-essence @) , the DE models including Chaplygin gas [7'j, agegraphic DE models [§, @], 
HDE model flol - [T3| and so on. 

The HDE model, which has got a lot of attentions in recent years, is based on Cohen et al. [l0| works who discussed 
that in quantum field theory a short distance cutoff could be related to a long distance cutoff (IR cutoff) due to the 
limit set by black hole formation. If the quantum zero-point energy density is due to a short distance cutoff, then the 
total energy in a region of size L should not exceed the mass of a black hole of the same size, namely L 3 pn < LM 2 . 
The largest L is the one saturating this inequality and therefore we obtain the HDE density as 0. 

p D = 3c 2 M 2 L~ 2 , (1) 

where c is a numerical constant and M 2 = 8irG is the reduced Planck mass. In the literatures, various IR cutoffs have 
been considered, which lead to various models of HDE [l5l-fi8|. Recently, Granda and Oliveros (GO) have proposed 
a new IR cutoff which includes time derivative of the Hubble parameter and is the formal generalization of the Ricci 
scalar curvature, L = ( aH + /Jff 2 )" 1 / 2 [19]. The advantages of GO cutoff is that the presence of the event horizon is 
not presumed in this model, so that the causality problem can be avoided [19j. Besides, the fine tuning problem can 
be solved in this model 0. 

On the other hand, the expansion rate of the Universe is explained by the Hubble parameter H = a/a, where 
a = a(t) is the scale factor of the Universe, while the rate of the acceleration or deceleration of the Universe expansion 
is described by the deceleration parameter, 
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However, the Hubble parameter H and the deceleration parameter q cannot discriminate various DE models since all 
DE models leads to H > 0 and d > 0 or q < 0. In addition, the remarkable increase in the accuracy of cosmological 
observational data during the last few years, compel us to advance beyond these two important quantities. A question 
then arise: how can we discriminate and classify various models of DE? In order to answer this question, Sahni, et. 
al., 0] and Alam, et ah, [5IJ, proposed a new geometrical diagnostic pair for DE. This diagnostic is constructed from 
scale factor a(t) and its derivatives up to the third order. The statefinder pair {r, s} is defined as 0 
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The statefinder pair is a geometrical diagnostic in the sense that it is constructed from a spacetime metric directly, 
and it is more universal than physical variables which depend on the properties of physical fields describing DE, 
because physical variables are, of course, model-dependent. Usually one can plot the trajectories in the r — s plane 
corresponding to different DE models to see the qualitatively different behaviors of them. For flat ACDM scenario the 
statefinder pair is {r, s} = {1, 0} [22:]. This fixed point provides a good way of establishing the distance of any given 
DE model from ACDM model. It was shown that the statefinder pair {r, s} can also be related to the EoS parameter 
of DE and its first time derivative |2C} ■ The statefinder pair is calculated for a number of existing models of DE having 
both constant and variable EoS parameter. It was argued that the statefinder can successfully differentiate between 
a wide variety of DE models including the cosmological constant, quintessence, the Chaplygin gas and interacting 
DE models (see [23|,[24| and references therein). In this paper we would like to study the statefinder pair parameters 
{r, s} for HDE in standard cosmology as well as DGP braneworld scenario. 

This paper is outlined as follows. In the next section we investigate statefinder for HDE in standard cosmology 
with considering Hubble radius and GO cutoff as systems’s IR cutoffs. We also plot the related figures which show 
the evolution of statefinder in each case. In section III, we extended our study to DGP braneworld and calculate 
statefinder pair {r, s} for HDE in a flat FRW Universe on the brane. The last section is devoted to conclusions and 
discussions. 


II. HDE IN STANDARD COSMOLOGY 


Consider a homogeneous and isotropic FRW which has the following metric 

dr 2 


ds 2 = — df 2 + a 2 (t) 
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where k = 0,1, —1 represent a flat, closed and open FRW Universe, respectively. Due to the fact that the density of 
baryonic matter is much smaller than the energy densities of DM and DE, we ignore the baryonic matter throughout 
the paper. Thus, Friedmann equation in spatially flat FRW Universe may be written as 


H = 


1 


3 M 2 


(,Pm + Pd)i 


where p m and po are the energy densities of DM and DE, respectively. From equation ([5]) we can write 

-f- Qd = 1. 

where we have used the following dimensionless energy densities definitions, 

Pm o PD 


— 


Hu = 


3M 2 H 2 ’ 3 M 2 H 2 ' 


(5) 


( 6 ) 


(7) 


Recent observational evidences provided by the galaxy clusters supports the interaction between DE and DM 25]. In 
this case, the energy densities of DE and DM no longer satisfy independent conservation laws. They obey instead 


Pm T 3 H Pm — Q, 

Pd + 3H(1 + u>d)pd = —Q, 


( 8 ) 

(9) 


where ojd = Pd/Pd is the EoS parameter of HDE and Q describes the interaction between DE and DM. The choice 
of the interaction between both components was to get a scaling solution to the coincidence problem such that the 
Universe approaches a stationary stage in which the ratio of DE and DM becomes a constant J2(|. Here, we choose 
Q = 3b 2 H(pD + p m ) as an interaction term with b 2 being a coupling constant. This expression for the interaction 
term was first introduced in the study of the suitable coupling between a quintessence scalar field and a pressureless 
cold DM field [Hill. 

First, we choose Hubble radius as IR cutoff, L = iJ _1 , thus the the hologeraphic energy density reads 


p D = 3c 2 M 2 H 2 . 


( 10 ) 


If we consider c as a constant, then from definition 0 the dimensionless DE density becomes a constant, namely 
£Id = c 2 . In this case, without interaction between DE and DM the accelerated expansion of the Universe cannot be 
achieved and we get a wrong EoS, namely that for dust, ujd = 0 [HJ. However, as soon as the interaction between 
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two dark component is taken into account, the accelerated expansion can be achieved and implies a constant ratio of 
the energy densities of both components, thus solving the coincidence problem [l3| . 

Taking the time derivative of the Friedmann equation (0 and using Eqs. ©-©, we get 


H 3/ ^ \ 

Jp = ^V 1 + w d^dJ- (11) 

Also, if we take derivative of Eq. (fTUl) with respect to time, after combining the result with (0 and ED, we obtain 
the EoS parameter as 
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The rate of deceleration/acceleration of the Universe expansion is characterized by the deceleration parameter that 
encoded the second derivative of the scale factor. In our case the deceleration parameter is found to be 


H 1 3 b 2 

H 2 ~ 2 ~~ 2(1 -c 2 )' 


(13) 


It is clear that in the absence of the interaction where b 2 = 0, we have wd = 0 and q = 1/2 > 0, which implies that we 
encounter a decelerated Universe. Furthermore, in order to find more sensitive discriminator of the expansion rate, we 
consider the statefinder parameters which contain the third derivative of the scale factor. The satefinder parameter 
r can also be expressed as 
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Taking the time derivative of both sides Eq. © and using Eqs. ©, ©-([HD, we can obtain the statefinder pair {r, s} 
as 


9 b 2 (c 2 + b 2 — 1) 

r = 1 H- 

2(1 -c 2 ) 2 

c 2 + b 2 — 1 
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The statefinder pair parameters {r, s} become constant depending on b 2 and c 2 parameters. For c 2 = 1 — b 2 we have 
{r, s} = {1,0}, which is exactly the result obtained for spatially flat ACDM model [22j. So we can find relationship 
between the coupling constant b 2 and HDE constant parameter c 2 . 

Next, we consider the GO cutoff as the system’s IR cutoff, namely L = ( aH + /377 2 ) -1 / 2 . This cutoff was first 
introduced by Granda and Oliveros (H} and can be considered as the formal generalization of the Ricci scalar curvature 
radius [3|. The cosmological implications of the HDE model with GO cutoff have been investigated in [29j. With 
this IR cutoff, the energy density © is written 

Pd = 3M 2 (aH 2 + (3H), (17) 

where a and /? are constants which should be constrained by observational data, and we have absorbed the constant c 2 
in a and /? in HDE density. For the case of GO cutoff, we do not need to take into account the interaction between DE 
and DM and without it we still have acceleration. Thus, for economical reason, here we consider the noninteracting 
case with b 2 = 0 = Q, so the continuity equations read to 


(18) 
(19) 

Now taking the time derivative of both sides of Eq. © and using Eq. ©, we get 


Pm T 3 H Pryi — 0, 

pu + 317(1 + lod)pd = 0. 


From Eq. ED one can obtain 
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FIG. 1: The evolution of fI_o versus redshift parameter z for HDE with GO cutoff in standard cosmology. Left panel corresponds 
to a < 1 and the right panel shows the case a > 1. 


The equation of motion for the dimensionless HDE density is obtained by taking the time derivative of relation 0 

( 22 ) 


n D = -2n D —. 

3M 2 H 2 H 


Combining Eqs. m and m with (E^l) . we arrive at 

n' D = (1 — fijy) 


3 + - a) 


(23) 


where a prime denotes derivative with respect to x = In a and we have used relation tto = H fl' D . The evolution of the 
dimensionless HDE density in terms of redshift parameter l + z = a^ 1 is shown in Fig. 0. From these figures 
we see that at the early time Qd —>• 0, while at the late time where a —» oo, the DE dominated, namely Qd —> 1. 
Combining Eqs. (ED, ® and (ED with (ED, we can obtain the EoS parameter as 


lod = — 


1 

flu 




The deceleration parameter can be expressed in this model as 


, H tt D -a 

q H 2 8 


(24) 


(25) 


When rip —> a = 1 we have u>d = —1 and q = —1, similar to the cosmological constant. Therefore, we investigate 
the behavior of wd for values close to a = 1. We consider the case with a < 1 and a > 1 separately. We find that for 
both cases, our Universe has a transition from deceleration to the acceleration phase around z « 0.6. The behavior 
of wd and q are plotted in figures 2 and 3. From these figures we see that for a < 1 the EoS parameter can cross the 
phantom line in the future, i.e., ujd < — 1, while for a > 1 the EoS parameter is always larger than —1. 

Taking the time derivative of Eq. (ED, after using (|23l) . we obtain 




(26) 


Inserting Eqs. ED and ED in ED, after some calculations, we obtain statefinder pair parameters {r, s} as follow, 
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FIG. 2: The evolution of u>u versus redshift parameter z for HDE with GO cutoff in standard cosmology. Left panel corresponds 
to a < 1 and the right panel shows the case a > 1. 




FIG. 3: The evolution of the deceleration parameter q against redshift parameter z for HDE with GO cutoff in standard 
cosmology. 


It can be easily seen that for a = 1, we have {r = 1, s = 0} as expected. The statefinder parameter r(z) is plotted in 
figure 4 and 5 for two cases. These results are compatible with the arguments given in 12(1 [2~ll . 


III. DGP BAREWORLD MODEL 

In this section, we would like to extend the study to braneworld scenario. According to this scenario, our Universe 
is realized as a 3-brane embedded in a five dimensional spacetime. A special version of braneworld scenario was 
proposed by Dvali- Gabadadze-Porrati (DGP) [3C}> in which our four-dimensional FRW Universe embedded in a 
five-dimensional Minkowski bulk with infinite size. The recovery of the usual gravitational laws in this picture is 
obtained by adding to the action of the brane an Einstein-Hilbert term computed with the brane intrinsic curvature. 
The presence of such a term in the action is generically induced by quantum corrections coming from the bulk gravity 
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FIG. 4: The evolution of the statefinder parameter r versus redshift parameter 2 for a < 1 (left panel) and a > 1 (right panel). 




FIG. 5: The evolution of the statefinder parameter r versus declaration parameter q for a < 1 (left panel) and a > 1 (right 
panel). 


and its coupling with matter li ving on the brane. It was argued that this term should be included in a large class of 
theories for self-consistency [31], |32|. In this model the cosmological evolution on the brane is described by an effective 
Friedmann equation which incorporate the non-trivial bulk effects onto the brane. The modified Friedmann equation 
in this model is given by [33] 



(29) 


where p = p m + pn is the total cosmic fluid energy 


density on the brane. The cross over length r c is given by 


Mpi = G 5 
2 Mf 2 G 4 ’ 


(30) 


and is defined as a scale which sets a length beyond which gravity starts to leak out into the bulk [33]]. In other 
words, it is the distance scale reflecting the competition between 4D and 5D effects of gravity. The parameter e = ±1 
represents the two branches of solutions of the DGP model [33]. For e = +1 branch, there is an accelerated expansion 
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at the late time without invoking any kind of DE or other components of energy [34j . however, it suffers the ghost 
instabilities problem [35]. Besides, it cannot realize phantom divide crossing by itself. On the other hand, for e = — 1 
branch cannot undergos an accelerated expansion phase without additional DE component }.37l ]. For iJ _1 <C r c (early 
time), the Friedmann equation in standard cosmology is recovered, 


H z 


3 M 2 ■ 


pi 


For the spatially flat DGP braneworld (k = 0), the Friedmann equation C51) reduces to 

e „ 1 


H 2 - H = 

ry 3 M„ 2 


(pm “b Pd)- 


We can rewrite this equation as 


where we have defined 


1 2e \J^lr c — “b 




1 


477 2 r 2 ' 

In the remaining part of this paper, we shall assume the continuity equations hold on the brane, 

Pm + 3 H Pm = 0, 

Pd + 317(1 + ujd)pd = 0. 


(31) 

(32) 

(33) 

(34) 


(35) 

(36) 


Following the previous section, we shall consider two cutoffs, namely, the Hubble radius and GO cutoff, as systems’s 
IR cutoff in the HDE density. First we consider L = if -1 . Inserting pc, = 3c! 2 M 2 H 2 into the Friedmann equation 
(O, we get 




(37) 

which can also be written as 

1 — c 2 — 2e^/fl rc = fl m . 

(38) 

Taking the time derivative of Eq. 

(l37l) after using Eqs. (l34l) and (l35l). we arrive at 



0 H _ -3H m 

H 2 1 — c 2 — ey / Q rc 

(39) 


We can also obtain the equation of motion for tt rc . For this purpose, we first take the time derivative of Eq. 
and then combining the result with Eqs. (1351) and (1551) . we arrive at 


From Eq. HMD we have 


_ — 3n r< .(l — c 2 — 2e v / H^~) 
1 - c 2 - e^/H rc 


Pd , 
wd = —-1. 


3 Hp D 

while taking the time derivative of pd = 3 c 2 M 2 H 2 , we find 

H 

Pd = 


(40) 


(41) 


(42) 


and hence the EoS parameter reads to 


e^/fl ra 


2 H 

WD ~~3 lP~ ~ 1 - c 2 - e^/TE^’ 


( 43 ) 










8 




FIG. 6: The evolution of wd and the deceleration parameter q versus redshift parameter z for the two branches of DGP 
braneworld. For both cases (e = 1, c < 1) and (e = —1, c > 1), the diagrams are quite similar to each other. 


where in the last step, we have used (l38l) and (I551) . For r c ^> 1 (f2 ro —>• 0), the effects of the extra dimension vanishes 
and hence we reach to the standard cosmology regime. In this case wjj —> 0, which is an expected result. This is 
due to the fact that in the absence of the interaction, the choice of L = U -1 leads to a wrong EoS, namely that of 
dust Q, and hence this choice for the IR cutoff in standard cosmology cannot leads to the accelerated expansion. 
However, from (1551) we see that in the DGP braneworld, the natural choice for the IR cutoff in flat Universe, namely 
the Hubble radius L = H~ 1 : can produce the accelerated expansion. This is one of the interesting results we find in 
this paper. 

The deceleration parameter is also obtained as 


H _ 1 36^0“ 

H 2 ~ 2 2(1-c 2 -ey 7 ^)' 


(44) 


The evolution of wjj and q are plotted numerically in figure 0 where we consider two branches of DGP braneworld 
e = 1 and e = — 1 with different values of c. In order to have reasonable behaviour for the evolution of the DE during 
the history of the Universe, we should take c < 1 in case of e = 1, while c > 1 in case of e = — 1. For both cases 
diagrams have similar behavior. From figure El we find that for e = ±1 we have a transition from deceleration to 
acceleration phase around z s=s 0.6. From figure El on e can see that EoS parameter wd gradually tend to —1 with the 
evolution of the Universe, while it is larger than —1 throughout the evolution of Universe, which is an indication for 
quintessence behavior [36| . 

Taking the time derivative of Eq. (1551) and using Eqs. (1551) and (1551) . we arrive at 


H ( H n 2 9(1 - c 2 - 2e v / H^J 

VtfV (1 -c 2 -eyn^j + 2(1-c 2 -ey/ShT)' 

Using Eqs. (l39l) . (1551) and (1551) . after some calculations, we obtain the statefinder pair parameters as 


_ , Q H i H _ 1 9e\/f2 ro (1 c 2 2e-^/f2 T . c ) 2 

r - +3 Jp + Jp- 1 4( 1 _ C 2 _ e ^J3 

r - 1 _ (1 - c 2 - 2eyTUj 2 

~ 3(g — 1/2) “ 2(1 -c 2 -e v / n^) 2 ' 


(46) 

(47) 


Next, we choose the GO cutoff for the HDE density in the framework of the DGP braneworld. We restrict our study 
to the current cosmological epoch, and hence we are not considering the contributions from matter and radiation by 
assuming that the DE dominates, thus the Friedman equation becomes simpler. Substituting 113 into Friedmann 
equation (1321) . we can obtain the differential equation for the Hubble parameter as 


H 2 {l~a)~— iJ-/3F = 0. 
r c 


( 48 ) 
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Solving this equation, the Hubble parameter is obtained as 


m = 


et J 

r c (l — a) + ciee^ 


( 49 ) 


where C\ is a constant of integration. Since for r c 1, the effects of the extra dimension should be disappeared and 
the result of jig], namely 


m = 


p i 


ot — It 

must be restored, thus the constant c\ should be chosen as [H] 

r c {a - 1) 


Cl = 


Substituting ci in Eq. AMD, we obtain 


H(t) = 


r c (l — a)(l — ) 


In term of the redshift parameter, we have 


H(z) = 


r c (l - a) L 


1 + (1 + zY 


(50) 


(51) 


(52) 


(53) 


Solving (1551) for the scale factor, yields 


a(t) = an ( — 1 


i{t) = a 0 ( 

Taking the time derivative of Eq. (|52j) . we arrive at 

e~P 


r■ 


H = 


H = 


Pr 2 (l — a)(l — ) 2 
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ee & r c 


| 2e' 3r c ^ 

F^v 1+ 1 _ e ih) 


3 2 r3(l-a)(l-e^<=)2 V ! - e e 


Therefore, it is easy to show that 
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In this case the EoS and the deceleration parameters are obtained as [34 


wd(z) = -1 - 


1 — a 


(1 + a)(l + z)-r- +2 


'-(l + (1 + z)~^){l + a(l + z)~^~)- 


(54) 

(55) 

(56) 

(57) 

(58) 

(59) 


q(z) = -1 - 


1 — a 


P 


1 + (! + «)' 


(60) 


To check the limit of (l59l) in standard cosmology where r c 1, we note from Eq. (l54l) and relation a/ag = (l + z) 1 
that (a/ao)^~ = (1 + z ) _ ? i —>• 0, as r c 1. Therefore, (1551) reduces to 


wd = —1 + 


2 a — 1 


3 p ’ 


(61) 
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FIG. 7: The evolution of wd versus redshift parameter z for HDE with GO cutoff in DGP braneworld. Left panel corresponds 
to a < 1 and f3 < 0, while the right panel shows the case a > 1 and /3 > 0. 




FIG. 8: The evolution of the deceleration parameter q versus redshift parameter s for HDE with GO cutoff in DGP braneworld. 
Left panel corresponds to a < 1 and /3 < 0, while the right panel shows the case a > 1 and /3 > 0. 


which is exactly the result obtained in [19]. It is important to note that the EoS parameter of HDE with GO cutoff 
in DGP braneworld scenario varies with time. Thus, it seems that the presence of the extra dimension brings rich 
physics. For a = 1 we have lod = —1, similar to the cosmological constant. The case with a < 1 and a > 1 should be 
considered separately. In the first case where a < 1, we find that the EoS parameter can explain the acceleration of 
the Universe provided (3 < 0. In the second case where a > 1, the accelerated expansion can be achieved for /? > 0. In 
both cases, our Universe has a transition from deceleration to the acceleration phase around 0.4 < z < 1, compatible 
with the recent observations f38| - l40) . and mimics the cosmological constant at the late time. In figures 0 and d8J) we 
have shown the evolution of the EoS and the deceleration parameter for HDE with GO cutoff in DGP braneworld. 

Combining Eqs. m and (15511 with © and m, after using Eqs. (15Hl and JUD: we obtain the statefinder pair 
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parameters as 


r = 1 + (1 — a) 


3j3(a~^~ + 1) + (1 - a) (2 + a~P~) 


/3 2 (a,—p- + l) 2 


s = 


1 — 1 (a — l) 2 


2(1 + o' 3 ) + 1 


6/3 


3/3 


(3/3(1 + o' 3 ) + 2(1 — q;)^ (1 + u 3 ) 


(62) 

(63) 


For a = 1 we have {r, s} = {1,0}. In term of the redshift parameter, 1 + z = a , with ao = 1 for present time, Eqs. 
and (l63l) takes the form 


r 


s 


1 + (1 

a — 1 
6/3 


a 


3/3((l + «)-?" + 1) + (1 - a)(2 + (1 + z)~P~) 


/3 2 ((1 + zy 


l ) 2 


(a- l) 2 _2(l + (l + z)V) + l _ 

3/3 (3/3(1 + (1 + 2 )^) + 2(1 - a))(l + (1 + z)^) 


(64) 

(65) 


The evolution of the statefinder pair parameters for HDE with GO cutoff in the framework of DGP braneworld have 
been shown in figures 151 1 121 Interestingly enough, from figure dill) we see that r diverge as q —> 1/2, which corresponds 
to the matter dominated Universe, and mimics the cosmological constant, namely r = 1 ,q = —1, in the far future 
where 2 —>• 0. 




FIG. 9: The evolution of the statefinder parameter r versus 1 + z for HDE with GO cutoff in DGP braneworld. 


IV. CONCLUSIONS AND DISCUSSION 

Following the discovery of the accelerated expansion of the Universe, many DE models have been introduced 
to explain such a cosmic acceleration. In order to discriminate various DE models, the statefinder diagnosis pair 
parameters {r, s} was introduced HU- This diagnostic is constructed from scale factor a(t) and its derivatives 
up to the third order. It was argued that trajectories in the r — s plane correspond to a wide variety of DE models. 
The spatially flat ACDM scenario corresponds to a fixed point {r, s} = {1,0} in the statefinder diagnosis diagram. 
Departure of a given DE model from this fixed point provides a good way for establishing the statefinder of the 
presented model from ACDM. 

In this paper, we have considered HDE in standard cosmology as well as DGP braneworld. We have chosen the 
Hubble radius L = H~ x and the GO cutoff L = (aH + j3H 2 )~ 1 / 2 , inspired by Ricci scalar curvature of flat FRW 
universe, as the system’s IR cutoff . We have investigated the evolution of the EoS, deceleration and statefinder 
parameters. We found that in standard cosmology for both proposed cutoff, the statefinder pair parameters {r, s} 
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z+1 


FIG. 10: The evolution of the statefinder parameter s versus 1 + z for HDE with GO cutoff in DGP braneworld. 



FIG. 11: The evolution of the statefinder parameter r versus deceleration parameter q for HDE with GO cutoff in DGP 
braneworld. 


become constant depending on the model parameters and with suitably choosing the parameters they restore those 
of ACDM model. However, in the framework of DGP braneworld, the situation become completely different and 
some interesting results can be observed. Interestingly enough, we found that although choosing L = H -1 as the 
IR cutoff, in the absence of interaction, leads to the wrong EoS of dust wd = 0 [131, here because of the effects 
of extra dimensions the natural choice for IR cutoff in flat Universe, namely the Hubble radius L = if -1 , can 
lead to an accelerated expanding Universe. In this case wjj and q , as well as the pair parameters {r, s} depend on 
Q rc = l/(4:H 2 r 2 ). For GO cutoff in DGP braneworld, the cosmological parameters can be obtained as the function 
of redshift parameter z, and in case of a = 0, we reach to {r, s} = {1, 0}, as expected. We have plotted the behaviour 
of {r, s} and observed that r diverges at the past as q —> 1/2, which corresponds to the matter dominated era, and 
mimics the cosmological constant with r = l,q= —1 in the future. 
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FIG. 12: The evolution of the statefinder parameter r versus s for for HDE with GO cutoff in DGP braneworld. 
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